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SUMMARY 

V. hen a rocket is launched, there Is a short initial period of ac- 
celeration during which the rocket is unstable. As the flight velocity 
increases, the aerodynamic forces acting on the fins and stabilisers 
become large enough to give stability. Various methods have been em- 
ployed to stabilize the rocket during this launching period. Guide rails, 
"zero length" launchers, booster rockets which produce high initial ac- 
celeration, and auto-pilot controlled nozzles are typical devices that 
have beon used. 

This is an investigation of the requirements of a nozzle control 
which would stabilize the rocket during the launching period. The con- 
figuration Investigated is unique in that the nozzle of the rocket is 
mounted as a compound pendulum, and the movement of the pendulum 
is utilized to furnish the signal for the uozzle control servo-mechan- 
ism, thereby eliminating tho need for gyroscopic elements in tho con- 
trol system. The pendulum motion of the nozzle caused by a change in 
flight attitude of the rocket is introduced into a computer which pro- 
duces an output signal proportional to the attitude of the rocket. This 
attitude signal is fed back to the nozzle control, which positions the 
nozzle . 

The results of the analysis showed that the rocket was unstable 
during the take-off period when the nozzle control acted on the rocket 
attitude signal alono. Stability over a narrow range of feedback gains 
was indicated for the system using a simple lead circuit as a nozzle 
control, or in other words, when the nozzle control acted on both the 
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attitude signal and the rate of change of attitude of the rocket. The 
damping characteristics of this system were poor. By changing the 
nozzle control function to include a response to the acceleration of the 
rocket attitude, the damping characteristics were improved and the 
range of feedback gains was widened. 
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I. INTRODUCTION 



The rockot stability problem has existed for many centuries. 
Historians record that rockets were first used in battle by the Chinese 
in 1232 A.D. 3y 1400, rockets were in widespread use as major wea- 
pon! by the warring nations of Europe. The poor stability of the rocket 
in flight was responsible for its gradual replacement by the cannon 
which gave more accurate trajectories. In the late 1700's, India em- 
ployed rockets with improvod stability characterises against the 
British. V hen British soldiers attempted to defeat the forces ofHydar 
/vli. Prince of Mysore, a trained corps of rockot bombardiers launched 
rocket missiles at the charging British cavalry, inflicting many casual - 
tics. The Indian rockets were made of iron cylinders about eight inches 
long and two inches in diameter, with a bamboo shaft about ten feet long 
attached to give stability. 

^.n English military officer. Colonel illiam Congreve, became 
interested in the Indian advances in rocket development, and worked 
to further improve the rocket for British use. Congreve devised a 
launching system consisting of a tube mounted at an anglo on a tripod 
base, similar to tho mortar. Long shafts were attached to the rockets 
to keep thorn flying in a straight lino. The British employed rockets 
against Napoleon at Boulogne, and in other bombardments with great 
success . 

During the ’ ar of 1012, Americans were first introduced to 
rocket warfare when the British bombarded Bladenburg, Maryland, and 
subsequently Ft. McHenry at Baltimore. An American by the name of 
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Wllliam Hale Improved the stability of the Congreve rocket by replac- 
ing the wooden shaft with three curved fins on the rear part of the 
rocket body. The exhaust gases were deflected by the fins, causing 
the rocket to spin in flight. These rockets were used in the war with 
Mexico, and were found to have superior trajectories to the stick sta- 
bilized types. Again the rocket declined in popularity partly because 
of its inferior trajectory compared to that of the gun. No other out- 
standing developments in rocket stabilization were made following 
Hale's work until early in this century. 

With the advent of the airplane, and an increase in knowledge in 
the field of aerodynamics, the stability problem of the rocket in flight 
became defined more clearly. The use of tail fins, or tail fins with a 
shroud ring as on the Bazooka, and later, the use of controlled airfoils 
solved the in-flight stability problem to a great extent. 

In the first part of the rocket's flight path there is a destabilizing 
effect from the change in mass of the rocket as fuel is burned. The 
aerodynamic forces must counter this tendency in order to stabilize 
the rocket. During the launching period, the aerodynamic forces are 
too small to counteract the effect of the change in mass, and the rocket 
is unstable or at best neutrally stable. Consequently, stability during 
the critical take-off period must be attained by some means other than 
the use of aerodynamic surfaces. 

Various means have been employed to stabilize the rocket during 
the launching poriod. Guide rails to hold the rocket to a straight flight 
path until it has gained sufficient speed to become aerodynamically 
stable have been used successfully. The length of the guide rail has 
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been reduced to a minimum in some applications by firing the rocket 
with a very high initial acceleration. The German V-2 rocket was 
stabilised during the launch by an automatic pilot which controlled 
four carbon vanes in the Jet stream behind the nozzle. The American 
"Viking" and the German "Ensian" employed a controllable nozzle 
positioned by an automatic pilot. 

A unique means of launching a large rocket with a booster rocket 
mounted free to swing as a compound pendulum was investigated by 
J. C. Norris in 1951 following a proposal by H. S. Tsien of Cal Tech. 
(Ref. 1). For the system investigated, the booster rocket was mounted 
below the main rocket on gimbals off the center of gravity of the 
boostor as shown in Fig. 1. A deviation of the main rocket from its 
course would cause motion of the booster relative to the main rocket 
causing a change in the direction of the thrust axis. It was hoped that 
a practical configuration could be devised to give launching stability. 
The results of the investigation were somewhat discouraging in that 
for a typical vehicle -booster combination for which stability was indi- 
cated, the required length of the tall boom was too great to be prac- 
tical. 

Since the basic idea of the pendulum mounted thrust unit appears 
sound, the possibility of utilizing the pendulum in conjunction with a 
servo-control is suggested. 

In this thesis, the requirements of such a servo-control system 
which would stabilize the rocket during the launching period are inves- 
tigated. The rocket has a controllable nozzle mounted as a compound 
pendulum; conceivably, the motor and nozzle of a liquid propellant 
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rocket could be mounted as a unit ae shown in Fig. 2. The motion of 
the pendulum with respect to the rocket is fed to a computer which cal- 
culates the attitude of the rocket. The attitude signal from the com- 
puter is fed to a nozzle control. The aim of this investigation is to 
determine the requirements of the nozzle control such that it will sta- 
bilize the rocket during the launching period. 

The equations of motion of the configuration were first derived, 
assuming all aerodynamic forces to be negligibly small during the 
launching period. Motion in one plane only was considered. External 
forces considered were the rocket thrust, assumed to be constant; the 
jet damping force arising from rotation of the jet stream; and the gra- 
vity force. The effect of the change in mass of the rocket resulting 
from the burning of fuel was considered. The displacements of the 
rocket and nozzle about the Initial condition were assumed to be small 
quantities for the period of the stability investigation. 

Following the derivation, the equations involving the pitching 
motion of the rocket and the motion of the nozzle were transformed to 
simple algebraic equations by use of the Laplace Transform. From 
the transformed equations, it was possible to solve for the transfer 
functions of all of the components of the control system, as shown in 
Fig. 4, except that of the nozzle control. The form of the nozzle con- 
trol transfer function was estimated, and the stability of the resulting 
system was investigated by use of the Nyquist Stability Criterion. The 
Evans Root Locus plot was used to find the damping coefficients of the 
modes of oscillation of the system. Successive estimates of the form 
of the nozzle control function were made in order to Improve the 
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atability of the ayatem. 

It waa found that for a nozzle control function employing feedback 
of the rocket attitude aignal only, the ayatem waa unatablc. For a con* 
trol function employing feedback of both tho rocket attitude and rate of 
change of attitude, it waa determined that a atable range exiated be* 
tween upper and lower limita of feedback gain; the limita depending on 
the conatanta of the control function. The damping of the ayatem waa 
poor . 

The nozzle control function waa reviaed by Introducing additional 
lead circuita, or aenaltivity to the accelerative change of attitude of 
the rocket. The damping waa improved by thia change, and the range 
of galna waa widened. 
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U. ANALYSIS 




Fig. 3. Block Diagram of a Feedback Control System 



A typical feedback control system which might be employed to 
guide a rocket of the controllable nozzle type is shown in Fig. 3 above. 
The attitude of the rocket, as indicated by the angle <jf } which the 
rocket makes with the vertical, is measured by an attitude gyro. The 
gyro signal is matched with the desired rocket attitude signal and the 
error, or difference, is fed to a nozzle control. The nozzle control 
applies a moment /t/ to the nozzle which results in a relative angular 
displacement (gy 0$ of the nozzle with respect to the rocket. The 
thrust axis rotates with the nozzle, causing a change in the rocket 
attitude . 

Now consider the case of a rocket with the motor and nozzle 
mounted as a compound pendulum as shown in Figs. 1 and 2. With any 
change in rockot attitude, the pendulum moves in a calculable manner. 
From a knowledge of the system dynamics and a time history of the 
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pendulum motion, the attitude of the rocket may be computed. Suppose, 
then, that the angular displacement of the nozzle from the center line 
of the rocket, y ^ , is measured and fed as an electrical or me- 

chanical signal into a computer with an output representing the attitude 
of the rocket . The gyroscopic element is thereby eliminated 
from the circuit. 
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Fig. 4. Block Diagram of a Feedback Control System 
Employing a Compound Pendulum Nozzle 

The block diagram of the proposed feedback control circuit as 

shown in Fig. 4 above may be interpreted as follows. 

/* *» 

The desired input of the coordinate y , or angle which the 
rocket makes with the vertical, is matched with a feedback signal rep- 
resenting the actual angle at any given time. The difference between 
the desired angle (in this case it is desired that ~ O for vertical 
takeoff) and the computed angle is fed into the nozzle control me chan- 
ism as an error signal, . The error signal initiates the 

/# ,i 

application of a moment ^ to the nozzle gimbals. The moment 
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contributcs to the nozzle dynamics to yield a certain position of the 
nozzle relative to the rocket, which is represented by ^ ^ . The 

thrust line is rotated off the center line of the rochet by movement of 
the nozzle; this contributes a moment and side force to the overall 
rocket dynamics, resulting in a final rocket position. 

The function! p z , p 3 , and aopearin in tho block diagram 
of JTi^. 4 are transfer functions derived from the equations of motion 
of the configuration which have boen transformed by the use of the 
L.& >lace Transformation. The transfer functions are linear operator- 
and can be defined by inspection of tho block dia ram of the system 

r Q« - Q3 tr _ c _ Qs 

(i . e . , \~i ~ , 13 — — Q 4 - Q-j )• The nozzle control trans- 

fer function p, is the unknown quantity which is the subject of this in- 
vestigation. 

The e juation relating the output of the system to the input may 
be derived as follows: 

Tollowin u the circuit from Q 3 )N to Q 3 0UT 




vround the loop 




F, F; F< Qj 1 * 
1 + f 



Substituting for Q-j in tho first equation 

A . . A nrn.iT Fl Fl 




F. F, F. j _ fl Ft F; 
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The System Transfer Function therefore becomes: 

K r — Output _ Fi_ ^ ^ 

~ Input [ + p, F* P 4 

The stability analysis thus involves 1) derivation of the equations 
of motion for the rocket and movable noazle configuration, 2) conver- 
sion of these equations by the Laplace Transform to solve for the trans- 
fer functions , and , and finally, 3) synthesis of a notzle 

control function P ( , and Investigation of the complete system trans- 
fer function G s f° r ^s stability characteristics. 

The Nyquist Stability Criterion and the Evans Root Locus Method 
are applied in the stability investigation. Refs. (2) and (3) contain ex- 
planations of the use of these and other methods of analysis applicable 
in the field of stability and control. 
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III. EQUATIONS OF l MOTION 



Tho assumed rocket configuration with all oxternal forces to be 
considered is shown In Fig. 5. A schematic diagram of the coordinate 
system and dimensions of the configuration are shown in Fig. 6. 

Frictional and aerodynamic forces are neglected. The change of 
mass of the main rocket body with burning of fuel is considered. Mo- 
tion in one plane only is considered. 

It is convenient to employ coordinates through in the 
derivation, but of these, only four coordinates, through <^ , are 

independent. The coordinates and ^measure the horizontal and 
vertical distances respectively to the center of gravity of the entire 
configuration, and have the dimensioni of length. Coordinates and 
measure the angular displacement of the main rocket center line 
and the nozzle line respectively from the vertical, and have the dimen- 
sions of the radian. 

The equations of motion in the four degrees of freedom, namely 
in the directions of coordinates \ ^ » and ^ arc derived 

using tho generalized Lagranglan momentum equation: 



sL (11 \ _ U = 2 

a<u~ 1 



(i) 



Elimination of Extraneous Coordinates 



Dy taking moments about the center of gravity of the system as 



shown schematically in Fig. 7, the relation between the system center 
of gravity and the centers of gravity of the main rocket body and the 
nozzle may be found: 
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Fig. 7. Sketch Showing the Relation Between the Centers 
of Mass, Mj and M 2 , and the Center of Gravity of the System 

z M c 3 . = H 3 (V ij -Mi 9 (vtJ= 0 

M, {%;\,) = M 2 (%,-%) <2 > 

From the geometry of the system as shown in Fig. 8 the follow- 



ing relations can be obtained: 
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Equating the vertical distances between mass centers: 



V% 5 ) 

0 - ^6 ~ % - 
M ~ co S ( V S) 

Equating the horizontal distances between mass centers: 

^5~ % - ^4 5 (%->, ~ <$) 

o - 

^ SLn( V <$) 

By considering similar triangles in Fig. 8 it follows that: 

<U - U - h ~ h 

V %> 



( 3 ) 



( 4 ) 
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M^-li) _ Ml 

Substituting Eq . (2) into (5), the result is 



= Mi(%-1a) < 6 > 

Eqs. (2), (3), (4), and (6) are four linear expressions containing 
the four extraneous coordinates ^ , %j> %g and the principle 

coordinates and . The former may be expressed 

in terms of the principle coordinates as follows: 

F rom Eq. (3) 



%(,- cos (t,, - 5) -t % Q 


(7) 


F rom Eq. (4) 




1-T = ' A S'- (^- S) 


(8) 


Substituting (8) into (2) 






(9) 


Substituting (7) into (6) 




*L.- \ 7 - M ' i 4 C05(«v.-5) 
re bt m.+m,. ^ > i / 


(10) 


Substituting (9) into (3) 






(ID 


Substituting (10) into (7) 






(12) 



It will be noted that the above equations contain the variables 



and $ from Fig. 7 which must be expressed in terms of the 
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princlpal coordinates. 

In Fig. 8 by the law of cosines 

\ZV~- t K (u) 

and by the law of sines 

j?, _ 4 j 

§ r 5 in" 1 

v i - k,- sin " l [|; 3 ,h CVt-b)] < u > 

To simplify Eqs. (13) and (14) consider the following approximations: 



CoS 



k~ I 



% t ~ Sin k 

Then from Eq. (13) 

JP 4 - t, - ft 

From Eq. (14) 

^3 ' ^ ' 17 ( % * ^3) 



^4 

q _ £ ~ q _ Al q 



and since 



4 = t -**. 
, - Jf. 



C '-3' % /-i, 



Tho approximations then are 
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Si "<v 

Cos C^"S) :r ' I 

A, - *i -**- 



(15) 



Substituting (15) into (9), (12), (11), and (10) yields the following re- 
lations between the extraneous and the principal coordinates: 



\ 5 - *V, +K-*M X 


(16) 




(17) 




(18) 


^8 = % t ~^= 


(19) 



The Kinetic Energy Term 

The kinetic energy appearing in the Lagrangian equation of 
motion may be derived as the sum of kinetic energies of translation 
and rotation by referring to Fig. 6: 

t = ±i,% +| l t % +i-‘) (20) 

“q '■ 

Derivation of the Equation of Motion for the Coordinate n 

The first term of the Lagrangian equation may be obtained by 
first taking the derivative of Eq. (20) with respect to ^ : 






H. 



+ i T 



Hi 
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The extraneous coordinates and appearing in the above expres- 
sion may be eliminated by differentiating Lqs. (16) and (18) and substi- 
tuting the results. Thus we obtain: 






Mi N, 






The above expression may be differentiated with respect to time to 
yield: 




The second term of the Lagrangian equation may be obtained by differ- 
entiating E{. (20) with respect to 



dr _ 

H, 



o 



(22) 



The generalized force 9 can be computed by considering the 
work done by the external forces as the coordinate Cj,. is varied, all 
other coordinates remaining fixed. 

In general. 



Work • — Si. ^ ^ 



or 




Work i 



In particular, for coordinate ^ , the force becomes, 

to I' ig. 6: 



3 



l 



Wor k | 




+ D CoSC^) ^9-/ 
4 S'. 



referring 



1 ith the approximations of E . (15), the above expression becomes 



a 






I 



F% + D 



(23) 
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Substituting Eqs. (21), (22), and (23) into Eq. (1), the equation 
of motion in the ^ ^ direction becomes, finally; 

* 1 - ^7 (*' ii - *'*) + 0$ (■*• V (l h) 



= H, + D 

Derivation of the Equation of Motion for the Coordinate 
In a similar manner, differentiating Eq. (20): 

+Miis h, 

and, differentiating Eq. (17): 

Hi - i 
1 

Similarly, differentiating Eq. (19) 

V 8 YZ - (M.-tMx) V J 

_ i 

H, 

Combining the above results: 

£1 

^ ' y M, + Mx V 7 






(24) 



By differentiating this expression with respect to time the first term 
of the Lagrangian equation becomes: 
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(25) 



(26) 



^ (|y = ( M , , MO r M, V; t u - ti ) 

The second term of the Lagrangian equation is: 

£-1 = o 

Hi 

The generalized force , referring to Fig. 5, may be ex- 
pressed as: 

~ _ Work* _ [-fM,4M x )9 + F Cos«^ - Dsin^] L 

«* 1 ~ — 

V ith the approximations of Eq. (15) 

2 1 c: •+ F - D (27) 

Substituting Eqs. (25), (26), and (27) into (1), the equation of motion 
in the ^ direction becomes: 



+ +M,i, F -3 (m,*m z ) - D<j^ 



(28) 



\ ' ^ •• 



Derivation of the Equation of Motion for the Coordinate ^.3 

The first term of the Lagrangian equation is obtained in an idem 
tical manner to that employed for coordinates ^ and : 
Differentiating Eq. (20): 



3 T _ . . 3 <i c 

» . — 1 1 7 + Mi — - 

M s 3 H: 






(29) 



and, differentiating Eq. (16): 



M, -+-fn t 
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also, by differentiating . (13): 



4 =. % - (i.P, )- MxM, (Q a - P. q ) 

>7 > < X> %/ ' °3 2 ^/ 



_ _ M. 

~ M , + M x 

Substituting into (29) 

IX 



f| f M,Ml Pi q 

|_ A * M,+Mj. J ^3 



M, Mi Qjz 



MitMi 






Differentiating the above expression \/ith respect to time, the first 
term of the L>agrangian equation becomes: 

♦ esrfrv**.! 



F rom Eq. (20) 



XL 

h % 



- "■(% 









6 5 












By differentiating Eq. (16): 









Mi. M, C, 

( M, + M») r 



y differentiating E*, . (18): 
^i 7 _ _ Mi. Mi -C i 

(M.+ M.) 1 



(30) 



(31) 



(32) 



(33) 



v i< • : n 

The term -1 which accounts for the propellent mass transfer out 

of the rocket is not correct, Inasmuch as this moment should be calcu- 
lated from the velocity with which the propellent leaves the transfer 
tubes and enters the motor. Corresponding errors appear in subse- 
quent equations. The numerical effect of the resultin terms is negli- 
gible . 
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In order to obtain ^ and 

sary to differentiate the exact expressions for ^ and , since the 
C^term has vanished in the linearization of the exact expressions. 
Eqs. (10) and (12) will be differentiated and the results then linearized. 
F rom Eq. (12) 

) ( e= %"^/ (12) 






±L 

^3 



it becomes neces- 



= % --£± - s<»e ® -Mi* t J2i- U cose 

^ MitKi 

^ * M» t-M^ M,tM v M.tMi o=y 3 

_ Mj_M. CoS 6 _ ^ 'dCQSQ J^liicose 

S^ 3 (/v\it Mi) S<y- 3 M.tM*. 'S ^ 



NU ^ coi© 



'^3 



Mi +Mt 

F rom Eq. (13): 

^4 - c<5 « C^V -S- 3 ) 

By differentiation: 

ai? 4 x ^ s in (%*-%*) 



^3 



■»< 



(30 



Mi- j, - -Ma^cvwJLii-t.] 

** " 4 a 

|4 = a j 0 ^t.„) 

^ 4 

£ & 



(35) 



(36) 



n- 



/ft / • • \ 

= - 1r (<v <-^) 



(37) 
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From Eq. ( 14 ): 

e = %- s '"~' [tz si ” (%-%)] 



By differentiation: 



e = % " 



( C V 4 -C kjj — ( 2 nd O •• t « r ■" ) 

y/~~\ ( Z nd ordtrT^r^vi) 



e - i 3 - -t (%- %) 



X 






30 

Bc b 

36 

Alto: 



.X 



( 38 ) 



( 39 ) 



X: 



%)] [- 



yi-[J±*:r(h’h)] £ ' 


" r 1 r a/ ^ ’ f T U/j | 


A± 


i ^ X 2 / — X i 

x " xT 




( 40 ) 


- CoS 0 -Ij- 

n 3 


^ X 1 


( 41 ) 


&C 0 S 9 __ s ,„ 8 |6 

a % 3< h 


X , 1 


( 42 ) 



By substituting Eqs. ( 35 ) through ( 42 ) into ( 34 ) and ignoring terms of 
second order or higher in the small quantities: 



H, 



a ^- M..M. '-3 

In a similar manner, from Eq. ( 10 ): 



Mi. 0 q , KiM, 0 q 



( 43 ) 



i = % + -^i 4 si« e e -£*- c«5e|| 4 - 



By differentiating this expression: 



-25- 



ils. z— [& *JUel5i!LP + linfl elk 

M.TMv L H 3 a ^3 H 3 



f CoSe|^ 31 M 0 — 

o<^ 3 a<y j 



M, [A 
|_H 5 



cos o- s ixie^A^. 

d *3 



+ 



!2ji. fi*l C038 -Sl.eP 4 l« 
i.*m,) x [3<v 3 Hj 



(M 



Again substituting Eqs. (35) through (42): 



(44) 



ir.M, • + 






He. = jn 



i<w. 



(45) 



Having evaluated and from the exact expressions, 

s t Hi 

0 can now be obtained by substituting these values in Eq. (31). 

Hi 



JlL 



Mi M, 



(M.tMj 5 

The generalised force 



K+m.. ^ 



(46) 



5, - 



Then, referring to Fig. 6 



- Work 3 

^3 



+ Ft,5 % ^[v^.^ 05 ^ +Dc °’%^r 3 [^r^‘ , ^ s, '“%J 



- D Si* 



U 
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Substituting in values of 9 and from Eqs. (11) and (10) 

7 6 

+ [F»s^- D *i»%] ^ [v 

Performing the differentiation indicated: 

S ^ lf - 1 , - 0-1.] [';$■/• 177 * S77 

Now, substituting Eqs . (35), (41), and (42): 

3i ~ ( 1 F V + 1 D ) {- fe. ^ + ( {- J] 

+ (F - 

Finally, retaining only first order terms, the expression reduces to: 

M . i, 



2 - - 
^3 M.t-M,. 



[F(^-h) + D] 



(47) 



The equation of motion for the coordinate can then be ob- 
tained by substituting Cqs. (30), (46), and (47) into Eq. (1). 

IA * M<) 

+ Ua 4, - H) 



(M.thq 3 >5 - KA [1^ W tu j 



(48) 
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V » M 

Derivation of the Equation of Motion for the Coordinate ^ 4 
By differentiating Eq. (20) with respect to: 

= M. + m, 4 |if + Mi 

■ rom Eq. (16): 

is - i, 

ifferentiating: 

ai-5 _ _ m..<i 

H* 

From Eq. (18): 



(49) 






Differentiating: 



Aii - 



k 



JPx 



f'di't M x 

abstituting these uantitios into Eq. (49) 

_ j ^ _ KJjJi 

1 r * M.*M» M T/ 

The first term of the Lagrange equation is obtained by differentiation 
of the above expression: 

*(!£) = J * % - (50 > 



F rom Eq. (20) 









■ y differentiating Eq. (16): 









Hi 



p 2 






(51) 
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0^5 — Ml M. $2- 



vnd by differentiating Eq. (18): 

M i 



i 7 = %. 



(M, 



_ Mi Mi 4 1. 
“ (M.tM,) 1 



In order to evaluate 



and 



* “Ve 



it is necessary, as 



in the case of the ^ coordinate, to differentiate the exact expressions, 
and then make simplifying approximations: 

F rom Eq . (12): 

cl - q _ Jll. \ + si m o 6 - Co5 e ^ £ 4 ccs6 



Differentiating: 




Mi jik.ee 



Mi jf o &lw Q 

M. ♦ Mj_ ^ 



e 






ji-e 



a© 

3<V^ 



_ Co5S 



M±iii a dCQ^o Ml ~°l±cose 
CHtMO -1 ^ H,tM L 



-v Mi ~ a c q 3 e 

Mi-tMu £ 



(52) 



By differentiating Eq. (13): 

^ ^ ^ J | ir - ^ 3 ) 
**4 ' JP* 



Also by differentiating Eq. (36): 

M± ~ hhii^hl 
hc )"> 
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11 y differentiating Eq. (14): 



de 

*>4 



ix 

^4 



And by differentiating E t . (38): 



3e _ -Mi (Vr'bX'k-'U') 




Then r<» may be evaluated by the substitution of the above ex- 

pressions and Eqs . (35) through (42) into Eq. (52): 



H© 

3 *4 



M i 

KtMj 




(t A.-iMO' > 4 



(53) 






F rom Eq. (10): 
M 






= C k 






i 1 * s'i* © © — 



M, ^4 

-i-4- coiflr;- 
M,tM t 2>-t 



- - 

(fA.tMv)' 



i»- 



cos e> 



Differentiating: 

a<v„ + t'iSlh* + s,'eeg± - 

k|_ H, H< ^ rt »Vt 

) La La r /N 4 






- cos© 



St 



?!_*)_ Mi Mi fn 4 

) (m.i-m») x I 3< W 



cose 



SO 

Hi j 



Then 



3c W 



may be evaluated in a similar manner to 



substituting for the extraneous terms in Eq. (54): 
d<Vfl _ „ M, l* • _ Mz. M.-l* 



Hi 



(54) 



by 



( 55 ) 



The second torm of the Lagrangian equation may be evaluated 
by substituting the derivatives obtained in the foregoing steps into 
Eq. (51): 
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dT M l A _ MiM-fi Q a 

Kt+rAv 





-J2i , . 

Y*‘ '3 


H,) + 


Ki 


M 


“Vi 


The generalized force 












a + - 


Wo r 


k 








*^4 






reference 


to Fig. 6, it is 


evident 


that: 








= <^4- D 


,cos ^0 


1 H, [ V 


-a+p- 


0 s %\ 




+• (Fcos^ - D 


,,MC h) 


T^, ['Vs ‘ 




) >,] 






- M,« 


1 

1 3<V4 







( 56 ) 



Substituting for the extraneous coordinates and simplifying, as in the 
derivation of Eq. (47): 




(57) 



Introducing Eqs. (50), (51), and (57) into Eq. (1), the equation of mo- 
tion for coordinate c p. becomes, finally: 

- •**h. -fifes*. . 



M. +Kj 



Mi M.ii 



l liliA (o ~o. \ q — - o (Ja^iSh 

( + “ \K+Mi 



(M 



-t- 






(58) 
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IV. COMPUTATION OF SYSTEM COMPONENT 
TRANSFER FUNCTIONS 



Only the two equations involving the rotation of the rocket (co- 
ordinate ^ ) and the motion of tho nozzle (coordinate ^ ^ ) are con- 

sidered since the stability investigation is made for rotation only. 
First, the thrust, jet damping, and moment of inertia terms are sub- 
stituted in terms of the equation variables and constants. 

Rocket thrust: 

p= — 3 I sp 

Jet damping force: 

D- - M,i s ^ 

(The negative sign is necessary in order to make the term 
positive, since M, is a negative quantity. ) 

Moment of Inertia derivative: 

I, - k'h *> radius of gyration) 

I,— M ,-k, 

• . r'l. ! 

X ( Cr K-^i (since will bo small for 

most configurations) 



Substituting the above terms in Eq. (48): 

M, Mi 






_ . fi _n \ ^ tr, C, \ - M.M.-fi • 



Mi Mi -$• 



M.-+M 



%^ 2 . M,+ Mi 



(59) 
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Again, substituting for p , D , and I t , and introducing the applied 
moment ^TlintoEq. (58): 




■f. K>i 

(M.i-NWp 



fih-J’-h) 



Mi M. -Pi 



(HrM.) 




MlMi 0 a 



+ £*i Mi J} L 

M.iMi 





+ th 

M.tMt 



i,) - *m. 



(60) 



The Laplace Transformation is applied to the equations prior to 

the stability analysis. Ref. (2) explains in detail the use of this trans' 

formation in such analyses. In essence, the transformation converts 

the differential equations to simple algebraic equations. 

• • 

Since the vertical and horizontal velocities and are 

• • 

quite small during the launch, the terms Involving ^ and ^ are 
neglected in the equations . 

Application of the Laplace Transformation to Eq. (59) yields: 






5 + 



_ M,tMi (M.tM ,) 2 



= o. 



+ [ frV M> ^ Mi M, $J± 

L MitMv J K\,-+ Mx __ 

j M.M.iglsf Mi 
L M.tMi (M* TmO 7 



(61) 



where represents the transformed coordinate , and S the 
complex variable. Similarly, transforming Eq. (60): 



The effect of the control moment "Vn, on the main rocket is neglected 
in Eq. (59). For a system in which the nozzle is of comparable size 
and weight to the main rocket, the term should be included. 
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l~\ J f M,M, 4. f» 1 c l [Ml M.-t.fl 1 q _ _ M J- 

Q 3 L'mttmt’ 1 s + L (m.+ sir J 5 - 1 V| (s> r 



0 . [Lv *$£ J ^ * [ 0 $- *Hu*&£)]s- [ 0 ^]j 

where |v) represents the transformed moment °ft\ . 

In order to solve for the transfer functions and it is 
necessary to rearrange the equations slightly; Eq. (61) becomes 

r\ r fl .* , f M m.M. f.JVI cl 

[r> A 1 JfKM.l.P.1 m wi'.rrl, .[nH.Sjhf 

[Wl yjJ[LTvTSTj S LO".’"O v KTSTJ 1 LH-iMv' (63) 

Similarly, Eq. (62) becomes: 

Qj - 1.] s ‘ + + *■ h ft *$&)] s ) 

- [Q,- o,] {[v ss£] s' * *■*&-£&)] s 

_fy£UjJ + ^ 0 

[.(m.+mo J j 



(64) 



Then the transfer function from Eq. (63) becomes 



F* = 



O 3 [ M.-tf'v J ^ T lSM.tMO x M.tMvJ (. MitMi J 

f -r i M,KliO t (g,-(vV j c' 1 ' ^ a + Mi M. £ -Pf 1 ! r 

L 1 K^Mi J L ' (Mi - ** 1 *.) 1 J 



(65) 



By substituting x-q. (63) into (64), the function may be written as 
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Q 4 - Q~s 

M 



[[ M.f 



< — 




* 

/ 



•is. 



l 




As an aid in the numerical computations for the stability analy- 
the transfer functions are symbolized as follows: 



where: 




Q 3 _ 

Q<^~ Q-ij 



p - Q<* - Qs _ 

“m 




(67) 



( 68 ) 



c V [ MfMA + , + fnhl.il,, . nlh'W l 

+ |^(Mr* Mi")' L M,+ Mv" (KtM 0 J Jj 









7 

I 



T 



Mt Mi v i ~ 0 ^ fv4,M, •(, fs 
(Mj tMi ) 1 M,iMi 





Hi O j I + | Mi &) 

Z J ^ L ( Mi +Mi) x 



+ 




1 * iWMi 

* [Cm,+ M 0 V 




NH ft. 
M.rMv 



)J S ’[ 



MLfvUll 



Wr 
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V. 3YNTI12TSI . TH1- NOZZLJi CONTROL. TR \NST£.R 




To Illustrate the determination of the nozzle control transfer 
function, a specific numerical example will be vorked out using the 
German V-2 as a typical vehicle. 

As indicated in the Analysis, the system transfer function may 
be written as 



as a numerical example, consider the V-2 rocket with pivoted 





(69) 



.•.Iso, the computor synthesizes the dynamics of the rocket to give 




(70) 



Substituting Lq. (70) into (69) 




(71) 



Let 



*6 



(72) 



.-ubstitutlng for p L and from Eqs. (67) and (68): 




(73) 



rocket motor and nozzle. The parameters X • Y • Z • anc ^ W 
have the following values as calculated in the .appendix: 
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X= ZQ7S S l - *6/S - /;o4f y ooo 



y= 72, 3 So s' - HO 5 



Z = Z o 3 o S z ~ //. 6 S S 



> 



( 74 ) 



fr/= 6 4-6 S z + /0.<5 > i‘S 

The problem now la to determine a function f~, which will give 
the system suitable stability characteristics during the launch. The 
procedure will be to: 

1. .(lume a function f-, . 

2. Check the system for stability by the Nyquist criterion. 

3. Having established a stable system, determine the roots of 
the system transfer function by means of the Evans Root 
Locus Method. 

4. Knowing the roots, at a typical operating gain of the feedback 
system, investigate tho damping of the system. 

5. If the resulting motion is unsatisfactory, improve the choice 
of Fi and repeat the process. 

In explanation of the methods to be used to investigate the sta- 
bility of the feedback control system, it will be noted that the equations 
of motion are linear differential equations with constant coefficients. 
The solution for the particular coordinate in question can then bo rep- 
resented in the general form 
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where 5 is a simple root of the equation. In order to have a damped 
motion which will subside with time, it is obvious that the complex 
roots s ( - must all have negative real parts, mother way of stating 
this requirement of a stable system is that all the roots s ( - of the sys- 
tem must appear in the left or negative half of the complex piano. 

The Nyquist Criterion 

The Nyquist Criterion is an application to tho simple feedback 

network of the Cauchy Theorem in the theory of complex variables. 

// •» 

This theorem may be employed to show that if is tho vector rep- 

resenting the function of a complex variable "s* , the function having 
m poles and / 7 " zeros within a closed contour C in the "s* plane, that 
as the point S moves around the contour once in a clockwise direction, 

// * // *1 

the vector C ^ carries out n-m clockwise revolutions about tho ori- 

// »» 

gin. If the contour C is chosen to include the entire right half piano, 
and the function in question is plotted as the complex variable encircles 
the right half plane, the difference between the number of zeros and 
poles of the function that exist in the right half plane will be indicated 
by the number of encirclements of the origin. Either the number of 
zeros, or poles existing within the contour must be known previously. 

In applying the theorem to the feedback network, the function in 
question is usually expressed as 

K <? = 

/ - rc,- 5j 

whore the function n can be factored; thereby determining the 

zeros of the function 4^ G s " . The denominator '/+ Sfc' in genoral 



cannot be factored so readily, but application of the Cauchy Theorem 
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will show the number of roots of the denominator existing in the con- 
tour chosen, or the right half plane. 

✓✓ * 

In order to obtain a plot independent of the gain K t the inverse 
of the function may be considered. 

-L = -L „( 

/r s 6 s /re 

/ / 

Then if the vector ^ q s encircles the origin, the vector 

y/ | I // 

will encircle the point -/ , as the variable s moves around the con- 
tour of the right half plane. Or similarly the vector -q~ will en- 

" j. “ " , •• .. 

circle the point 'T . The function -£~ is then plotted as s en- 
circles the right half plane. By application of this Nyquist Criterion, 
the stability or Instability of the system can readily be determined. 

yv \\ 

If the nozzle control transfer function is first assumed to be 
equal to a constant ( F - / - K t ) the system transfer function takes the 
form 



fsG s = 



* '6 



/ + ATG 






xz - yw 



7 



(75) 



For the Nyquist diagram, the inverse of G is plotted as the vec- 
tor 5 encircles the right half of the S plane. The function to be 
plotted then is 



/ _ xz- y*v_ 

G (1) X 



(76) 



In order to use the Nyquist stability criterion it is necessary to know 
whether the function A 'G has any zeros in the right half plane. By 
factoring Eq. (73) the following expression is obtained: 
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_ # r „ n -dl ( Asv '9.79] fs- /?.&// 

*g, = A/- nix /o yJ (77) 

' ( S fs-.006jfs-(-. O2,0i+t6.2sj£s-(- ,oZ0*/-Ci>.2SjJ 

Since the function Ac''has one aero In the positive real half of the 
plane, the Nyquist Criterion requires that the plot of q~~ must 

encircle the point - once in a counter-clockwise direction as the 
S vector circles the right half of the s plane once in a clockwise 
direction. 

Fig. 9 is the Nyquist diagram of Eq. (76). .vs the frequency 
increases from zero to positive Infinity, thence to negative Infinity and 
back to zero, there is one clockwise encirclement of all gains from 
O to oo . The net encirclement may be considered to be t.»o clock- 
wise rotations, indicatin. two roots of tho system transfer function 
existing in the right half plane. The system is unstable for all gains. 

Next, assume the nozzle control transfer function to be of the 
form £= 0+ / S) * then 



Ac 6, = 



f i f/r 7s) {"y * -y/v) 






(78) 



The Ny juist diagram for the above equation is shown in Figs. 10 and 
11. xgain the number of zeros of existing in the right half plane 

determines tho Nyquist criterion. Tho factored equation now become t> 



-clff s - {f)X T + '« 79 ^ s - /? - 6/ J 



(V ) 



*fysr rr tM* 2 >*M n r 

y-S- ood s-fr02oy, t&zsMsf-azoy- Stir)! 

Since there is one zero in the right half plane, again it is desired to 

M M 

have one counter-clockwise encirclement of the gain o as an indica- 
tion of stability. 
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It was found necessary to plot the Nyquist diagram to several 
scales to indicate clearly the crossing points of the negative real axis. 
The enlarged portions were plotted for the case when 7 =/ , and ap- 

pear as Figs. 12 and 13. From these plots it is clear that there is 
one counter-clockwise encirclement of all gains from /T*- /-2 to 
Af' - 7$~ . The configuration is stable for gains in this range. 

The Evans Root Locus Method 

This method is a graphical solution for the variation of the roots 

'' » *•" 

of the feedback system transfer function G as the gain K is 
varied from aero to infinity, where 









/ * p fsi 

The plot is constructed by observing the following rules (Ref. 2): 

" II 

1. For n z O , the poles of 4T S G s are the poles of 

2. For A'-—- c$ , the poles of ^ q" are the zeros of 
with any missing poles taken at s = cq • 



f * ,| 

3. The product of the distances from a root to the poles oi n g 
divided by the product of the distances from the root to the 

A 1 1 

zeros of A' G is equal to the gain 4 ' •' 

,, 

4. The sum of the angles from the zeros of /f " 6 to a root, less 

** ii 

the sum of the angles from the poles of 7T& to the root must 

/ * ** A It 

oqual 7 ? , or 7? plus an integer multiplied by ^ 77 , 



- ii 

The locus of the roots of /f' £ can be determined graphically 

5 <5 

by trial and error, observin. the above rules. The gain scale can be 
plotted along the path of the roots, and for a chosen operating gain. 
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thc root* can be determined by inspection. 

The root locus plot of Eq. (18) for the case when 7=-/ ie shown 
in Fig. 14. The roots emerge from the poles of ^ G at aero gain and 
proceed along the loci Indicated to the scros of ^ C at a gain of infin- 
ity. Gains at intermediate points can bo approximated by dividing the 
product of the distances to the poles of by the product of the dis- 

tances to the scros, and by the 'constant gain factor in the equation 

(49.2 x 10 -6 ). It is apparent from tho root locus plot that the range 

* * " 

for a stable configuration is from the gain where the roots s 

// »» 

cross into the stable half plane to the gain where the roots y 
cross into the unstable half plane, referring to Fig. 14. The loci in 
this range are plotted to an enlarged scale and appear in Figs. 15 and 
16. As an example of the computation, approximate expressions for 
the gain along the loci in this range are: 

For root 




~ .. .~a>. , .. 1 ^ < .. .. 



For root 

,3 CL £ ^ (for ~J small) 

C' 1 ,/ 5T77')f> / sy. 7* T&sYi/J?. u 39/)y?.z x/o* 

3 £70 fri- 02 0</) 

It will be noted that the gains at the limits of the stable range 
check within limits of accuracy of the scale used on the graph as shown 
in the table below for the Nyquist and the Evans Root Locus plots. 
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Method 


Upper Limit of K 


Lower Limit of K 


Nyquist plot 


75.0 


12.0 


Evans Root Locus 


74.80 


11.85 



Having established the gain scale along the root loci in the range 

of positive stability, referring to Figs. 15 and 16, it is possible to 

^ * 

choose an operating gain and determine the values of the roots ^3 

* 

and y at that gain. Some possible operating points have been 
tabulated below 



Roots 


K = 50 


K 3 60 


K « 70 


Si 


-.0031 + i .156 


-.00325 + i . 171 


-.0033 + i .185 


S 3 


-.0031 - i .156 


-.00325 - l . 171 


-.0033 - i . 185 


s 2 


-.0068 + 1 6.25 


-.004 + i 6.25 


-.0003 + i 6.25 


s 4 


-.0068 - i 6.25 


-.004 - i 6.25 


-.0003 - i 6.25 



" vN 

The low frequency roots Sy ^ are associated with the motion 

of the main rocket, whereas the high frequency roots Sy y are as- 
sociated with the motion of the nozzle. The time to damp to half am- 
plitude and the period may be computed for the modes quite simply; 
the results of such a computation for a gain of A'* £ O are tabulated 



below. 
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II 

o 

o 


Rocket Oscillation 


Nozzle Oscillation 


Time to damp to 
half amplitude 


210 seconds 


171 seconds 


- eriod of one 
cycle 


36. 3 seconds 


1 second 


Number of cycles 
to half amplitude 


5 . 7 cycles 


171 cycles 



It la apparent that the damping la quite poor, and an improve- 
ment ia desirable. The acceleration of the V-2 at take-off ia roughly 
equal to 30 ft. /sec.^. If it is assumed that at a velocity of 200 ft. /sec 
the aerodynamic forces stabilize the rocket, the period of take-off in- 
stability ia approximately 7 seconds long. For the foregoing example, 
the rocket conceivably could complete only one -fifth of a cycle of 
slightly damped motion in the take-off period as a result of a unit dis- 
turbance at the instant of firing. The nozzle oscillations are so rapid 
as to conceivably cause overloading of the control device. 

In order to improve tho dampin.. characteristics of the system, 
referring to the root locuo plot. Fig. 14, the addition of one or moro 
zeros along the negative real axis would serve to pull the roots S, 
and Sj into the negative half plane more quickly, since zeros attract 
the roots much as lines of force are drawn to magnetic poles. This 
should increase the magnitude of the negative roal part of the roots at 
a given gain. Also, the angle of departure of the root would be 
increased counter-clockwise as a result of the angle condition of the 
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root locus plot. The root would then be forced to travel a longer path 
before passing into the unstable half plane. 

The addition of teros in the left half plane involves adding addi- 
tional lead networks, which in real components have transfer functions 

of the type C/r£JS) where € » To 

Assuming a new nozele control transfer function to be of the form 



/r _ 7~s j f/+ rs)f /+ rs) 

5 J 0\ +jf 0 Sj 



where 7~= / , tho factored expression for /C Q is now 

) (83) 

(S(S-00 <oJ(S /2oJfSr 2 0)k-(- C20V,c ' 6 . 2^jfs-{zOZ CV-c '6-25jj J 

The root locus plot of Eq. (S3) is shown in Fig. 17. The range 
of gain for stability is now from K = 4 to K ® 918. 

Dy approximating the loci of the roots in the stable range by 
straight lines, the roots associated with various gains can be estimated 
conservatively. Some representative values are tabulated below. 



Root 


K = 600 


Si 


-.0446 + i .541 


s 3 


-.0446 - i .541 


S 2 


-.02 + i 15.25 


s 4 


-.02-1 15.25 



If the damping time and period are computed for the latter 
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system, the following results are obtained. 



K = 600 



Rocket Oscillation 



Nozzle Oscillation 



Time to damp to 
half amplitude 



15.3 seconds 



34. 2 seconds 



Period of one 
cycle 



11.6 seconds 



. 41 2 seconds 



Number of cycles 
to half amplitude 



1 . 3 cycles 



33 cycles 



Interpreting the results as before, the rocket now completes ap- 
proximately two-thirds of a cycle of oscillation durin the launchin 
period, but the effect of a unit disturbance is damped to a much greater 
extent. The nozzle oscillations occur at a higher frequency but th** 
damping has been improved considerably. 



VI. CONCLUDING REMARKS 



The analysis shows the feasabillty of stabilising a rocket during 
the launching period by utilising the motion of a compound pendulum 
suspended nozsle as the input to a feedback servo-control system. It 
was determined that the nozzle control for such a system must include 
response to both the attitude and the rate of change of attitude of the 
rocket, as computed from the pendulum motion, in order to obtain 
stability. It was further determined that changing the nozzle control 
to Include response to accelerative motions of the rocket improved the 
damping characteristics of the system. 
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APPENDIX 



Tho V-2 Rocket 






Wt. in Kg. 


Wt. In Lbs. 


Moment Arm, ft. 


Moment 
lb. ft. 


1. Explosive 


930 


2160 


7.42 


16, 000 


2. Fuselage 


1750 


3860 


26.65 


103, 000 


3. Pump 


450 


993 


37.5 


37,200 


4. Motor 


550 


1214 


41.7 


50,600 


5. Auxiliaries 


300 


662 


14.0 


9, 300 


6. Alcohol 


3500 


7720 


20.83 


161,000 


7. Oxygen 


5250 


11600 


30.2 


250, 000 


Totals 


12, 980 


28, 620 




727, 100 



Moment arm of C.G. = ^ ^ = 25.4 ft. 

28, 620 

Note: Data was obtained from Ref. (4) and is only approximately 
correct. 
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Calculation of tho Moment of Inertia 




Anume the rocket is a solid cylinder as shown In the above sketch. 
Radius of gyration of a solid cylinder about its C.C. is 

a /h iLlEJ ' 

3 r v L 



Additional Data from Ref. (4) 

F = Total Thrust =27.2 tons = 54, 400 lbs. 



W = Total Weight = 12. 980 Kg. = 28, 620 lbs. 



Mi * Mass Flow =- . U i 320 lb8t » - 193.2 lb. /sec. = - 6 lb. sec. /ft. 

100 sec. 



F 

= -i 

Ml E 



54, 400 
193.2 



= 282 see. 



I 8 p = Specific Impulse 



5 92 
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Calculations for the Constants of Section 111 



M, = 2 7, 406 I bs = 651 



I b 5 ec 



ff 



M, - 1,2 1 4 I b 5 = 37.7 

_ m-? 5 lb - l l _ 5ec 

M| - 193-2 5ec - ^ j- 

kr. - 9.7 ft 

|^ t ae 2. ft ( r 4 I u S of jyritiorv of motor) 

i, = 20 ft 
7*. 3.7 ft 

7, * 0.3 ft 
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